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Abstract. Heavy-quark production provides a sensitive probe of the gluon structure of nucleons 
and its modication in nuclei. It is also a key probe of the hot-dense matter created in heavy-ion 
collisions. We will discuss the physics issues involved, as seen in quarkonia and open heavy-quark 
CNI ■ production, starting with those observed in proton-proton collisions. Then cold nuclear matter ef- 

fects on heavy-quark production including shadowing, gluon saturation, energy loss and absorption 
will be reviewed in the context of recent proton-nucleus and deuteron-nucleus measurements. Next 



we survey the most recent measurements of open-charm and J / y/s in heavy-ion collisions at RHIC 
and their interpretation. We discuss the high-pr suppression and flow of open charm in terms of 
. energy loss and thermalization and, for J /yr, contrast explanations in terms of screening in a decon- 
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fined medium vs. recombination models. 

Keywords: quarkonia, gluon saturation, quark gluon plasma 
PACS: 24.85,+p, 25.75.-q, 25.75.Nq 



CO 

O : CHARM PRODUCTION IN P+P COLLISIONS AT RHIC 
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duced state and how it hadronizes remain uncertain. Absolute cross sections can be re- 



produced by NRQCD models that involve a color octet state[l], but these models predict 
transverse polarization of the J /yr at large pj that is not seen in the data|12j]. A general 
complication in understanding J/y/ results is the fact that ~40% of the 7/ y/s come from 
decays of higher mass resonances {yr' and #c)L3-] - a feature that may contribute to the 
lack of polarization seen. One exception to this feature is the maximal transverse polar- 
ization observed for the states[4]; where the lack of feed-down for these states 
may allow the polarization to persist. 

J/y/ cross section measurements for p+p collisions at i/j = 200 GeV from 
PHENIX||5|] are shown in Fig. [TJ These results, based on approximately 500 J /y/s 
from the 2003 run, provide the baseline for both CNM studies in d+Au collisions and 
QGP studies in A+A collisions at RHIC, and are presently one of the limiting factors 
in obtaining precise nuclear modifications. However p+p data from the 2005 and 2006 
runs will soon improve this baseline significantly with over 40,000 J /y/s. 

Open charm measurements at RHIC suffer from large systematics and statistical un- 
certainties due to the statistical subtraction methods that are used. Measurements by 
PHENIX and by STAR differ substantially [(J on the size of the charm cross section. 
Also the measured cross sections lie substantially higher than current theoretical predic- 
tions as shown in Fig. El 
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FIGURE 1. J/\j/ cross section vs rapidity 
for 200 GeV p+p collisions at RHIQ5]. 
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FIGURE 2. Open charm plus beauty 
cross section from prompt electrons di- 
vided by FONLL theory vs vr for 200 
GeV p+p collisions at RHIcfl]. 



NUCLEAR EFFECTS ON CHARM 



When quarkonia are produced in nuclei their yields per nucleon-nucleon collision are 
known to be significantly modified. This modification, shown vs. xp in Fig. |3] for 800 
GeV p+A fixed target measurements and in Fig.H|at RHIC energy, is thought to be due 
to several CNM effects including gluon shadowing, initial-state gluon energy loss and 
multiple scattering, and absorption (or dissociation) of the cc in the final-state before it 
can form a// y. 
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FIGURE 3. Nuclear moification fac- 
tor a vs xp for and y/' pro- 
duction in ^fs = 38 GeV collisions 
in E866/NuSea|7], and for D° from 
E789H. 
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FIGURE 4. Rapidity dependence of the 
J/y/ nuclear modification factor, R^Au for 
200 GeV d+Au collisions at RHIC If]. 



Shadowing is the depletion of low-momentum partons (gluons in this case) in a 
nucleon embedded in a nucleus compared to their population in a free nucleon. The 
strength of the depletion differs between numerous models by up to a factor of three. 
Some models are based on phenomenological fits to deep-inelastic scattering and Drell- 
Yan dat alllOl). w hile others obtain shadowing from coherence effects in the nuclear 
mediumlllll Il2ll . In addition, models such as the Color Glass Condensate (CGC)lTj] 



yield shadowing through gluon saturation pictures where the large gluon populations at 
very small x in a nucleus generate a deficit of gluons at small x. 

In the final state, the produced cc can be disassociated or absorbed on either the 
nucleus itself, or on light co-moving partons produced when the projectile proton or 
deuteron enters the nucleus. The latter is probably only important in nucleus-nucleus 
collisions as the number of co-movers created in a p+A or d+A collisions is small. 
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FIGURE 5. Test of scaling vs X2 andxF for J / \j/ suppres- 
sion data for three different collision energies. Data is from 
Refs.0E@ 
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FIGURE 6. Nuclear dependence of 
heavy quark suppression vs pj from sin- 
gle muons in PHENIX. 



However, J/\jf suppression in p(d)+A collisions remains a puzzle given that one does 
not find a universal suppression vs X2 as would be expected from shadowing, Fig. |5k; 
while vs. Xf the dependence is similar for all energies, Fig. Od. This apparent xp scaling 
supports explanations that involve initial-state energy loss or Sudakov suppression lflUl . 

On the open-charm front, there are no substantial modifications seen at central rapidity 
in d+Au collisions, but for forward rapidity (shadowing region) - as shown in Fig. In- 
substantial suppression is seen, while some enhancement is see at backward rapdity 
(Au-going direction). 



J /iff IN HEAVY-ION COLLISIONS - A QUARK GLUON PLASMA 

SIGNATURE? 

One of the leading predictions for the hot-dense matter created in high-energy heavy- 
ion collisions was that if a deconfined state of quarks and gluons is created, i.e. a quark- 
gluon plasma (QGP), the heavy-quark bound states would be screened by the deconfined 
colored medium and destroyed before they could be formed llal . This screening would 
depend on the particular heavy-quark state, with the l// and Xc being dissolved first; 
next the J/y and then the T's only at the highest QGP temperatures. The CERN 
SPS measurements!!^ showed a suppression for the J/y and 1// beyond what was 
expected from CNM effects - as represented by a simple absorption model constrained 
to p+A data. In addition to explanations involving creation of a QGP, a few theoretical 
models B 1 VII were also able to explain the data without including a QGP, so the evidence 
that a QGP was formed was controversial. 



The first measurements from PHENIX at RHIC in 2004 are beginning to yield results 
- see Fig. [5] for preliminary results for Au+Au and Cu+Cu collisions 11 911 . First it is 
important to understand what the normal CNM J/xy suppression should look like in 
these A+A collisions. This is illustrated by the blue error bands for A+A collisions 
in Fig. [U which represent identical theoretical calculations to the analogous blue error 
band in Fig.|7]for d+Au collisions. As can be seen the present d+Au data lacks enough 
precision to provide a good constraint on the CNM effects. As a result it is difficult to 
be very quantitative about the amount of "anomalous" suppression observed in A+A 
collisions, although there does seem to be a clear suppression beyond CNM for the most 
central collisions. 
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FIGURE 7. Results for J[w sup- 
pression in d+Au collisions|5] com- 
pared to a theoretical calculation 
that includes absorption and EKS 
shadowing l20ll . 




FIGURE 8. J ' jx\t suppression in Au+Au and Cu+Cu 
collisions for forward rapdidity and central rapdityl 19] 
compared to predictions for CNM from the same cal- 
culations as shown in Fig. 020]. 



On the other hand, all of the models JTvL EH l22ll that were successful in describing the 
lower energy SPS data over-predict the suppression compared to the preliminary data at 
RHIC - unless a "regeneration" mechanism is added as was done by Rapp ll22ll and by 
Thews[23]. The regeneration models assert that if the total production of charm is high 
enough then densities in the final state will be sufficient to have substantial formation 
of 7/ty/s from the large number of independent charm quarks created in the collision. 
This production mechanism was almost insignificant at SPS energies but at RHIC may 
be substantial. This leads to a scenario in which strong screening or dissociation by a 
very high-density gluon density occurs to a level of suppression stronger than the RHIC 
data shows, but the regeneration mechanism compensates for this and brings the net 
suppression back up to where the data lies. This is shown in Fig.|9l 

An alternative interpretation of the preliminary results, sequential screening, is given 
by Karsch, Kharzeev and Satz 124fl . In this picture, they assume that the 7/i// is never 
screened, as supported by recent Lattice QCD calculations for the 7/ \y - not at SPS nor 
at RHIC. Then the observed suppression comes from screening of the higher-mass states 
alone (y/7 and Xc) that, by their decay, normally provide ~40% of the observed 7/i//s. 
This scenario is consistent with the apparently identical suppression patterns seen at the 
SPS and RHIC shown in Fig. HO] 

As a result we are left for the moment with two different scenarios that provide 
explanations for the RHIC A+A data. Both include the QGP in their picture, either 



through color screening in the QGP or through severe suppression of the J /\\f by a very 
high gluon density. Further tests from the data will be necessary to clarify the picture. 
Regeneration models predict narrowing of both the rapidity and pj distributions, but so 
far the preliminary data shows little or no change in the rapidity shape from ordinary p+p 
and only a hint of narrowing of the pr- We are also trying to extract a measurement of 
flow for the J / ty/, since emerging results for single charm are beginning to show flow and 
the 7/va's, if they were from regeneration, would inherit this flow. These tests await the 
more precise final analysis of the 2004,5 Au+Au and Cu+Cu data; and higher statistics 
runs for Au+Au and d+Au in the near future. 
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FIGURE 9. Theories that agree with SPS data FIGURE 10. Universal dependence on energy 

do not agree with RHIC data, unless regenera- density of J/y/ suppression measurements at 

tion is added as in the Rapp (solid blue) curve. RHIC and at the SPS[24]. 



OPEN CHARM IN AU+AU COLLISIONS 

In Au+Au collisions, open charm (and beauty) together, are suppressed due to energy 
loss in the dense medium, with gluon densities per unit rapidity of up to 1000 infered 
in some theoretical analysis. However, as shown in Fig. HU calculations that include 
both radiative and collisional energy loss ll25ll predict too small a suppression when both 
charm and beauty are included. Flow has also been oberved for heavy quark production. 
As shown in Fig. [T2J the flow is similar to that of light quarks at small pr, but at 
higher pj the data with large uncertainties hints at vanishing flow, consistent with 
simple expectations that higher pj charm simply punches out of the medium and never 
thermalizes. 

SUMMARY 

Substantial uncertainties remain in the understanding of charm production cross sec- 
tions, and the polarization of charmonia. There are also a number of cold nuclear matter 
effects that influence their production in nuclei and cloud our understanding of the sup- 
pression seen in nucleus-nucleus collisions. Two competing pictures are able to explain 
the J/\]f suppression seen in nucleus-nucleus collisions at RHIC - one involving se- 
quential screening in the plasma of the various charmonia states; the other with strong 



dissociation of all charmonia states by a dense gluon field but recombination of inde- 
pendently produced charm quarks. For open charm, the the energy loss observed in the 
dense medium from nucleus-nucleus collisions is larger than that expected from theoret- 
ical models that include radiative and collisional energy loss of both charm and beauty. 
Higher statistics data with higher luminosity runs as well as RHIC vertex detector up- 
grades will enable more precise data in the future that will give a clearer understanding 
of the rich physics in charm production. 
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FIGURE 11. Energy loss calculations 
compared to open heavy (charm + beauty) 
data vs pt !25Ti . 



FIGURE 12. Elliptic flow of open heavy 
(charm + beauty) compared to Rapp 
calculations|261. 
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